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The  Application  of  coupled- mode 
theory  to  propagation  in  shallow  water 
with  randomly  varying  sound  speed 

Rolf  Thiele 


Abstract:  Sound  propagation  in  the  Baltic  Sea  cannot  be  described  cor¬ 
rectly  by  range  independent  propagation  models,  even  for  areas  with  more 
or  teas  constant  water  depths.  The  reason  is  the  inabihty  of  inch  models  to 
take  acconat  of  the  effect  of  forward  scattering  of  sonnd  produced  by  horl- 
sontai  variations  in  the  speed  of  the  sound.  Schneider  and  Selbchopp  have 
simalated  the  forward-scattering  effect  by  applying  a  parabolic  equation 
method  in  conjunction  with  a  Monte  Carlo  ray  tracing  program  (MOCAS¬ 
SIN).  This  paper  treats  forward  scattering  in  the  normal  mode  concept, 
using  the  mode  coupling  formnlation.  given  by  S.T  McDankf.  The  mode 
coupling  was  added  to  the  SUPERSNAP  normal-mode  code.  Comparisons 
are  made  with  the  Schneider-Sellachopp  calculations  and  with  real  data. 

Keywords:  Baltic  ,  S'  conpled-mode  theory  e  MOCASSIN  o  'mode 
coupling  modelling  a  Monte  Carlo  ray  tracing  program,  normal  mode 
theory  S'  propagation'  Scattering  st  shallow  water  'o~STJaP  o> sound 
speed  rf  source  localisation  methods  A  transmission  loss ,  ^ 
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1.  Introduction 

There  are  several  basic  methods  which  can  be  used  to  solve  the  wave  equation  for 
sound  propagation  in  the  ocean.  However  for  shallow  water  it  is  the  normal  mode 
technique  that  is  generally  preferred  [1].  This  means  solving  the  partial  differential 
equation  by  separation  of  the  variables  in  an  orthogonal  coordinate  system,  thereby 
converting  to  a  set  of  ordinary  differential  equations.  The  boundary  value  equation 
in  the  vertical  direction  forms  an  eigenvalue  problem.  Each  eigenfunction  is  an 
independent  solution  representing  an  individual  wave  propagating  with  a  specific 
phase  velocity  and  attenuation. 

The  basic  theoretical  assumption  for  the  use  of  the  normal  mode  technique  is  the 
horizontal  invariance  of  the  ocean.  This  rigorous  requirement  in  the  theory  may 
be  weakened  to  slow  horizontal  variations  through  the  adiabatic  approximation, 
which  still  neglects  the  transfer  of  energy  from  one  mode  to  another  (i.e.  the  mode 
coupling).  Milder  [2]  and  Pierce  [3]  reformulated  the  separation  for  slow  horizontal 
variations  with  local  eigenfunctions  and  included  the  effect  of  mode  coupling. 

McDaniel  [4-6]  used  the  Pierce  formulation  to  calculi  te  the  coupling  of  inodes  by 
forward  scattering.  Her  formulation  was  directed  to  horizontal  variations  of  the 
sound  speed.  Nevertheless  she  reduced  the  application  to  horizontal  variations  of 
the  boundaries. 

Schneider  [7,8]  has  shown  that  under  certain  conditions  the  forward  scattering  ari¬ 
sing  from  variations  in  the  sound  speed  may  have  a  strong  influence  on  the  sound 
transmission  loss.  Therefore  for  the  Baltic,  where  sound  propagation  is  very  sen¬ 
sitive  to  forward  scattering,  it  is  worthwhile  using  the  McDaniel  mode-coupling 
formulation. 

The  change  of  the  transmission  loss  observed  in  the  Baltic  Sea  is  used  to  verify 
the  mode  coupling.  However,  even  where  the  change  of  transmission  loss  is  not 
significant,  mode  coupling  may  influence  strongly  the  coherence  of  the  sound  field. 

in  principle  the  given  method  replaces  the  calculation  of  the  sound  field  for  range- 
dependent  sound-speed  profiles  (which  are  not  known  deterministically)  by  a  per¬ 
turbation  method  that  uses  a  range-independent  modal  computation  for  the  mean 
sound-speed  profile.  The  statistical  expectation  of  the  horizontal  gradients  of  the 
sound  speed  profile  is  used  for  the  calculation  of  the  expected  mode  coupling. 
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2.  Theory 

The  wave  equation  for  harmonic  time  dependence  is 

+  *2V  =  0,  (1) 

where  V  is  the  wave  velocity  potential  and  the  wave  number  k  =  2t//c  may  vary 
arbitrarily  with  depth  due  to  variation  in  the  speed  of  sound  e.  The  speed  of  sound 
is  assumed  to  be  slowly  varying  in  the  horizontal.  We  separate 

N 

=  51  *n(*>y)“n(* !  *.y)  (2) 

n=  1 


with  ♦„  independent  of  z,  where  un  is  a  local  eigenfunction  of  the  z-separated  part 
of  the  wave  equation 


dz* 


+  ( k 1  -  kl)  u„  =  0. 


(3) 


The  function  u„  can  be  chosen  such  that  the  orthogonality  relation 


for  n  =  m; 
otherwise 


(4) 


is  fulfilled.  Then  substituting  Eq.  (2)  into  Eq.  (1),  and  multiplying  by  and 
integrating  over  z  yields  (2): 


d'*„ 

dz*  dy* 


+  flmn 


9*m  r 

a  i  *'mn  o 
oy  OX 


}• 


(5) 


where 


Amn  -  dxi  +  gy}  J  dz , 

t" 

Bmn  =2  J  U„ 


dun 
'  dy 


Cmn=2  j\m^-dz. 


(6) 

(?) 

(8) 


are  coupling  coefficients  dependent  on  z  and  y. 
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We  propose  for  simplicity  a  propagation  direction  for  x,  such  that  di/dy  =  0. 
Further,  d2u„/dx!  is  assumed  small  compared  to  diin/9xi.  This  means  Amn  =  0 
and  Bm„  —  0,  reducing  Eq.  (5)  to 


a»»n 

dz 1 


+  *£*„ 


m^n 


(9) 


For  computational  reasons  McDaniel  [5|  assumes  a  gaussian  correlation  function  of 
the  coupling  coefficient 


(<7mn(*  +  A*)-Cmn(*)>  =  (Cln)e-**''L‘,  (10) 

where  Lx  is  the  coupling  correlation  length.  This  assumption  is  discussed  later. 
The  intensity  of  a  single  mode  (1)  is 


Pn(z,z) : 


WP  1  ttnU.KU)  flfa+jo.,. 


4c  V  2x z  y/k^ 


(11) 


The  coupling  may  be  described  as  the  change  in  power  of  each  mode,  given  by 


p 

dx  2-  v"m  ™ 


(12) 


with 


Pn  =  /  Pn(x,z)dz. 

Jo 


Perturbation  theory  gives  [4j  for  small  coupling  and  n  /  m, 


Power  conservation  may  be  taken  into  account  by 


Qnn  —  ^  ^  Q mn  —  2Qfn, 

mjtn 


(13) 


(14) 


where  a„  are  the  attenuation  coefficients  of  the  modes  without  coupling  (i.e.  volume 
attenuation  in  the  water  and  in  the  bottom). 

Furthermore,  from  McDaniel  [5]: 


C  - 

i-  mn  — 


kl 


3(*a) 

’  dz 


u.  dz. 


(15) 
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To  simplify  the  integral  it  is  assumed  k„  km  s:  k  and  c  a:  Co ,  so  that 


C  mn  — 


2k 


(*m  -  *n)fO 


dc  A 
,^u„dr 


(15a) 


may  be  used. 


The  variance  of  the  horizontal  gradient  of  the  sound-speed  grad^c  =  dc/dr  is  the 
basic  oceanographic  input  for  realistic  calculations.  This  has  to  be  depth  dependent, 
as  otherwise  the  orthogonality  integral  would  yield 


dc 

dz 


dc  fH  , 

fjx  J  ® 


for  n  /  m, 


with  the  result  that  coupling  would  be  neglected. 


In  reality  grad^c  should  have  a  rather  short  correlation  length  in  z.  When  this  cor¬ 
relation  length  is  small  compared  to  the  'vertical  wavelength'  of  the  mode  functions 
the  double  integral 


(f  rtin( *  4-  Ar)C  mn(^)) 


2k 


(In  )C() 


n,  .  ,  ,  /0e(*  +  A*,z)dc(*,r  + r')\  .  , 

um{z)um(z  +  i  )  ^ . ~dt - ^ - J  u„(z)u„(r  +  z)dzdz 

may  be  reduced  to 

<o  = 


2k 


(bn  bm)co 


1  J  ,H 

Lz  Um(2)un(2K8radlc(i))di-  (16) 

JO 


Eq.  ( 16)  can  be  applied  directly  to  Eq.  (13),  giving  the  coefficients  that  are  constant 
in  r  for  the  linear  system  of  differential  equations  of  the  first  order  Eq.  (12). 


This  system  of  differential  equations  forms  a  standard  real  symmetric  matrix  ei- 
genproblem 

(Q-AE]P  =  0,  (17) 

which  is  solvable  by  available  library  routines. 


The  resulting  N  eigenvectors  Pn,  each  with  a  constant  cm  and  relative  contribution 
coefficients  cm„  of  the  modes  m  with  c’m  =  1,  form  a  set  of  independent 

solutions  for  the  modes 


N 
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The  coefficients  c,„  are  determined  by  inversion  of  the  matrix  of  eigenfunctions,  and 
using  the  initial  values  of  the  mode  distribution 

UJO  1  -> 

Pmo  =  a  (19) 

4c  2*k,n 


For  large  distances  the  eigenvector  with  the  smallest  attenuation  will  remain,  gi¬ 
ving  the  stationary  solution  with  precisely  this  attenuation  and  range-independent 
vertical  field  distribution  of  the  associated  eigenfunction:  The  relation  of  the 
eigenvalues — which  are  attenuation  coefficients — may  be  used  to  determine  a  range 
at  which  the  cylindrical  spreading  law  additionally  to  the  mentioned  attenuation 
will  be  applicable.  Therefore  this  calculation  will  be  a  contribution  to  the  everla¬ 
sting  discussion  of  whether  shallow  water  propagation  is  best  described  by  a  10  log  r 
or  a  151ogr  range  dependence. 


3.  The  correlation  function  of  the 
horizontal  sound-speed  gradient 

The  most  comprehensive  description  of  the  variability  of  the  internal  structure  of 
the  ocean  is  the  internal  wave  model  of  Garrett  and  Munlc  (9).  This  model  gives  a 
spectral  slope  of  -2.5. 

Other  authors  |10,11|  propose  a  wavenumber  decay  with  an  exponent  of  -2  (instead 
of  -2.5)  as  more  likely.  Horizontal  differentiation  effects  a  multiplication  of  the 
horizontal  wavenumber  spectrum  by  the  factor  it2,  resulting  in  spectra  for  which 
the  Fourier  transforms  are  not  performable  bv  Rietnann  integrals  [13] 

This  means  that  McDaniel's  assumption  of  horizontal  correlation  in  Eq.  (10)  is  not 
justified,  and  consequently  we  cannot  be  sure  whether  the  frequency  dependence  of 
the  calculation  is  true 

Also,  as  this  is  a  first  approach  there  has  to  be  a  discussion  of  how  to  choose  the 
horizontal  correlation  distance  Lx,  when  we  know  that  it  may  in  reality  be  zero 
The  effect  of  L,  on  the  coupling  coefficients  is  twofold:  it  yields  a  scaling  factor  of 
the  variance  of  the  horizontal  gradient  and  in  the  exponential  part  of  Eq.  (13)  it 
yields  a  spatial  low-pass  filter  for  the  modal  interference  This  second  effect  reduces 
the  coupling  of  those  modes  with  the  larger  differences  of  order;  in  other  words  it 
limits  the  angular  spread  of  the  associated  rays  Since  the  coherence  length  should 
be  shorter  than  the  distance  of  modal  interference  (between  a  few  meters  and  some 
100  m)  it  is  appropriate  to  eliminate  this  filtering  effect 
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The  degeneration  of  a  parameter  with  its  horizontal  coherence  length  may  indicate 
a  change  in  the  power  law  of  the  associated  process,  i.e.  a  change  in  the  wave¬ 
number  Equation  (13)  exhibits  a  frequency-independent  solution,  which  applies 
for  large  scale  scatterers  (12],  for  small-scale  scatterers,  a  second-order  dependence 
(Rayleigh  scattering!  is  applicable 

1'ntil  now  no  actual  measurements  of  the  horizontal  gradients  were  available  for 
comparison  The  assumed  correlation  lengths  where  justified  by  comparison  of  the 
mode  conversion  calculation  using  Baltic  Sea  propagation-loss  experimental  data 
with  the  parabolic-equation  (PE)  method  of  Schneider  [8|. 


4.  Coherence  loss  by  mode  conversion  I 

Each  mode  has  a  certain  phase  velocity  and  all  modes  have  a  common  initial  phase. 

When  the  environment  is  known  exactly,  the  interference  pattern  of  the  modes  can 
be  determined  exactly,  and  conversely,  with  known  modal  structure  the  source 
position  can  be  determined 

This  deterministic  phase  relation  is  lost  by  the  mode  coupling  process  for  two 
reasons 

•  The  difference  in  the  phase  velocities  of  the  modes  causes  the  phase  of  a 
converted  mode  to  depend  on  the  horizontal  position  of  the  conversion  area 

•  In  the  vertical  direction  the  modes  change  in  phase  by  180°  at  each  node, 
which  means  that  the  positions  of  the  nodes  of  different  modes  are  different 
and  independent  Consequently,  the  resulting  phase  depends  on  the  depth 
at  which  the  coupling  takes  place 

One  consequence  of  the  mode  conversion  process  is  therefore  a  loss  of  phase  infor 
illation  the  converted  energy  is  incoherent,  while  the  unconverted  part  of  the  sound 
field  remains  coherent. 

The  calculation  of  mode  conversion  can  therefore  yield  an  evaluation  of  coherence 

loss.  This  is  valuable,  because  to  estimate  the  coherence  loss  by  directly  measuring  | 

variations  in  the  sound-speed  profile  is  not  easy 
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5.  Numerical  examples 

Two  examples  are  given  here  of  the  use  of  the  mode  conversion  (MC)  process: 
case  A  gives  the  initial  fit  in  comparison  to  Schneider  and  Sellschopp  (8j  and  case  B 
gives  a  typical  example  with  a  sound-speed  profile  that  represents  the  Baltic  duct. 
In  both  cases  the  real  variations  in  the  horizontal  gradients  were  unknown. 


Case  A.  Figures  1  to  4  show  the  rough  range-depth  distributions  to  be  in  good 
agreement.  However  the  computations  with  and  without  horizontal  variability  show 
significant  differences  in  the  shapes  of  their  contour  lines.  The  reasons  for  these 
differences  are: 

•  In  the  PE  process  the  modes  are  added  coherently  and  the  interference  has 
to  be  smoothed  out  by  averaging  in  range  and  depth,  whereas  in  the  MC 
process  the  modes  are  added  incoherently  without  averaging  in  range,  and 
only  in  certain  situations  is  there  some  averaging  in  depth. 

•  In  the  MC  process  the  assumed  stationary  variations  of  the  sound-speed 
profiles  are  used  to  derive  a  mean  solution,  whereas  in  the  PE  process  an 
experimentally-determined,  individual  two-dimensional  sound-speed  profile 
is  used. 

•  For  these  rather  high  frequencies  Schneider  was  forced  to  run  the  PE  process 
with  the  minimum  acceptable  resolution  in  order  to  save  computer  time. 
Thus,  some  numerical  inaccuracies  are  to  be  expected. 

•  Possibly  there  is  some  error  in  the  MC  process  as  a  consequence  of  the 
conversion  to  the  continuous  part  of  the  modal  spectrum  being  neglected. 

In  addition  the  computation  without  coupling  and  mode  conversion  shows  signifi¬ 
cant  differences  in  the  absolute  values  of  the  transmission  loss;  these  can  not  be 
explained  at  present.  A  calibration  error  in  the  PE  calculation  is  a  possible  reason. 

The  type  of  depiction  used  is  rather  sensitive  to  errors,  since  the  volume  attenuation 
is  neglected  in  both  calculations.  As  a  result  the  residual  attenuation  is  low  and 
errors  are  relatively  enhanced. 

What  is  remarkable  is  that  the  energy  scattered  into  the  duct  is  markedly  frequency- 
independent,  and  is  close  to  100  dB  at  about  25  km.  This  is  in  accordance  with 
Mellen  et  al.  [14],  who  suggest  that  the  loss  attributable  to  forward  scattering  will 
be  frequency  independent .  The  energy  in  the  duct  is  only  attenuated  by  scattering 
and  geometrical  spreading,  because  volume  attenuation  is  not  introduced  and  no 
direct  bottom  interaction  takes  place, 


Case  B.  This  case  deals  with  the  comparison  of  the  MC  process  with  reality, 
whereas  case  A  concerned  a  comparison  with  another  computational  method  using 
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the  same  data. 

Relative  to  case  A,  the  assumed  horizontal  gradients  had  to  be  increased  by  a 
factor  of  10  before  good  agreement  was  achieved.  This  is  surprising,  since  the 
temperature  profiles  of  case  A  were  measured  in  a  similar  area  of  the  Baltic  as  the 
acoustic  measurements  of  case  B 

Three  explanations  for  this  difference  are  offered: 

1 1 )  The  absolute  calibration  error  referred  to  above  leads  to  a  wrong  lit.  Increa¬ 
sing  the  average  horizontal  gradients  in  case  A  would  increase  the  scattered 
energy  by  approximately  10  dB  This  is  still  within  the  calibration  difference 
of  the  PE  process  and  the  inode  conversion  process  without  coupling. 

(2)  The  horizontal  spacing  of  the  measurements  for  the  temperature  profiles  was 
about  100  m.  In  the  PE  process  a  linear  interpolation  is  introduced  between 
the  measured  profiles.  To  get  correct  results  the  spacing  of  the  profiles  has 
to  be  shorter  than  the  horizontal  correlation  length  of  the  temperature  field 
or  possibly  of  the  horizontal  gradients  of  it.  Unfortunatly  the  correlation 
length  of  the  horizontal  gradients  is  expected  to  be  shorter,  and  then  the 
interpolation  means  a  low-pass  filtering  of  the  original  field  which  reduces 
the  average  horizontal  gradients 

(3)  Knowledge  about  the  changes  of  the  horizontal  gradient  field  in  complicated 
areas  such  as  the  Baltic  is  inadequate.  In  this  environment  gradients  will 
be  produced  more  by  changes  of  the  salinity  and  temperature  of  the  water 
masses  in  a  density  balance  than  by  internal  waves  Therefore  changes  in  the 
average  gradients  by  a  factor  of  10  are  not  unlikely,  since  the  measurements 
are  not  taken  at  the  same  position 

Figure  5  shows  the  comparison  with  the  measured  data  The  shallow  receivers  exhi¬ 
bit  the  effect  of  the  mode  coupling  significantly.  The  range- independent  solution  is 
completely  wrong  for  this  situation:  with  increasing  frequency  a  difference  in  the 
stationary  attenuation  arises,  presumably  because  the  volume  attenuation  in  the 
low  salinity  water  was  not  taken  correctly  into  account  In  the  model  computation 
the  attenuation  is  reduced  proportional  to  the  salinity  according  to  the  Thorpe 
attenuation  |15|,  but  the  correct  low  salinity  attenuation  is  still  controversial  |16) 
The  data  show  the  parallelism  of  the  transmission  loss  in  the  duct  and  beyond  as 
being  in  agreement  with  the  mode  coupling  The  sparing  between  these  parallel 
curves  is  correct  for  all  frequencies 

In  rase  A  it  was  shown  that  the  amount  of  scattered  energy  is  almost  frequency 
dependent  But  this  is  not  the  case  for  the  change  in  the  vertical  structure  of 
the  transmission  loss  In  Fig  6a  at  250  Hz.  (here  is  nearly  no  influence  by  mode 
coupling  at  10  km,  and  at  50  km  only  the  difference  outside  the  duct  is  significant 
At  500  Hz  (Fig.  6b)  the  fine  structure  in  the  duct  is  smoothed  out  at  10  km,  while 
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the  general  tendency  of  the  vertical  transmission  toss  does  not  change  even  at  10  km. 
At  1  kHz  (Fig.  6c)  this  change  in  the  fine  structure  is  observed  at  3  km;  at  10  km 
the  structure  changes  significantly. 

With  increasing  frequency  the  transmission  loss  profile  is  smoothed  by  the  coupling, 
giving  a  range  independent  structure  (Figs.  6d  and  e)  which  is  a  projection  of  the 
sound-speed  profile  (Fig.  13).  This  may  be  observed  best  at  the  small  change 
in  sound  speed  at  65  m  and  the  strong  vertical  gradient  at  90  m.  Between  the 
corner  points  of  the  sound-speed  polygon  which  are  depicted  in  the  transmission 
loss  profile,  the  transmission  loss  changes  very  smoothly. 

The  reason  why  the  energy  distribution  in  the  scattering  process  seems  to  be  almost 
independent  of  frequency  while  the  vertical  structure  of  the  transmission  loss  is 
strongly  dependent  on  frequency,  becomes  evident  when  one  looks  at  the  mode 
spectra.  Figures  7  to  10  show  the  energy  distribution  of  the  modes,  each  with  their 
ray  angle  corresponding  to  the  minimum  value  of  the  sound  speed. 

In  Figs.  7  to  10  the  range  of  1  m  is  only  given  as  a  reference  of  the  mode  distribution 
for  the  source.  The  values  were  extrapolated  from  the  far  field,  {is  the  far-field 
solution  of  the  Hankel  functions  is  also  applied  for  this  short  range.  The  figures 
were  obtained  via  the  plotting  program  for  the  mode  intensity  versus  angle  of  the 
SNAP  program  (1).  However,  it  is  the  complete  mode  energy  in  the  given  range 
that  is  depicted,  and  not  the  intensity  of  the  modes  at  a  single  receiver  depth. 

Figure  7b  (the  uncoupled  result)  reveals  some  relation  to  the  distribution  of  Fig.  7a. 
There  are  significant  differences  between  the  individual  modes,  and  the  high-order 
modes  are  strongly  attenuated.  The  coupled  modes  have  a  rather  smooth  envelope. 

Figure  8a  shows  the  initial  distribution  with  the  source  outside  the  duct.  In  Fig.  8b 
it  is  obvious  that  the  peculiarities  of  the  initial  distribution  are  better  retained  by 
the  uncoupled  solution;  however,  some  of  the  initial  structure  is  visible  in  both  the 
coupled  and  the  uncoupled  results. 

Figures  9a  and  b  show  that  the  source  depth  has  a  much  stronger  influence  at 
higher  frequencies.  For  the  uncoupled  situation  the  source  pattern  is  maintained 
for  all  modes  up  to  mode  40,  i.e.  the  propagation  effects  a  ‘low-pass  mode-filtering' 
with  a  cut-off  at  approximately  the  40th  mode  (Fig.  9c). 

With  the  receiver  outside  the  duct  only  the  modes  above  the  100th  are  excited 
(Fig.  9b).  Therefore  only  residual  modes  are  left,  and  the  mode  spectrum  is  thinned 
out  (Fig.  9d). 

With  mode  coupling  (Figs.  9e  and  f)  we  find  a  completely  different  situation.  The 
envelope  of  the  mode  spectrum  is  rather  smooth,  with  nothing  left  of  the  initial 
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pattern.  The  figures  that  have  different  source  depths  are  only  distinguishable  by 
the  different  common  levels  of  all  modes. 

Figures  10a  to  d  show  the  results  for  different  frequencies  and  source  depths.  The 
similarity  of  the  envelopes  is  evident.  The  envelope  reflects  some  of  the  relations 
between  the  sound-speed  profiles  and  the  attenuation  in  the  bottom.  Its  main 
feature  is  similar  for  all  frequencies  but  it  is  contorted  with  frequency. 

The  paradox  of  the  level  of  scattered  energy  (case  A)  seeming  to  be  frequency 
independent  while  the  structural  changes  are  strongly  frequency  dependent  may 
now  be  answerable.  The  level  of  scattered  energy  is  determined  by  the  envelope 
of  the  mode  spectrum.  The  envelope  is  only  weakly  frequency  dependent  for  the 
coupled  mode  result  but  strongly  so  for  the  uncoupled  result.  The  higher  modes 
are  heavily  attenuated  at  high  frequencies.  Equation  (16)  reveals  that  the  coupling 
between  the  modes  is  proportional  to  the  squared  reciprocal  of  the  relative  wave 
number  difference  corresponding  to  the  absolute  angular  difference.  Therefore  the 
coupling  is  strong  at  high  frequencies,  with  the  angles  close  together,  and  the  fine 
structure  of  the  mode  spectrum  is  smoother  at  higher  frequencies. 

Figure  11  shows  examples  of  the  transmisson  loss  over  range  and  depth  similar  to 
those  for  Figs.  1  to  4  (case  A).  There  is  a  remarkable  similarity  between  Figs,  lid 
and  e,  which  relates  to  different  source  depths.  The  sound  speed  profile  is  again 
duplicated  in  the  contour  lines. 

The  coherence  of  the  sound  field  is  demonstrated  by  Fig.  12.  Although  the  vertical 
structure  of  the  coherence  depends  on  the  source  depth,  the  range  of  each  individual 
contour  line  is  on  average  independent  of  the  source  depth.  However,  the  coherence 
can  be  said  to  depend  strongly  on  the  frequency,  since  the  neighbouring  modes 
interact  strongly  at  high  frequencies.  This  does  not  affect  the  transmission  loss 
significantly,  but  the  coherence  is  lost  whenever  there  is  mode  conversion. 
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6.  Conclusions 

Forward  scattering  attributable  to  sound-speed  variations  strongly  affects  propaga¬ 
tion  in  the  Baltic.  This  was  shown  by  Schneider  and  Sellschopp  using  the  parabolic 
equation  (PE)  process  [8]  and  is  here  demonstrated  by  a  mode  coupling  (MO) 
method. 

Variable  sound-speed  profiles  have  a  strong  effect  on  both  the  transmission  loss 
and  the  coherence  of  propagation.  The  forward  scattering  can  not  only  increase 
the  attenuation  but  also — with  the  source  or  receiver  beyond  the  optimum  depth— 
reduce  it. 

This  report  only  treats  forward  scattering  caused  by  sound-speed  variations.  Essen 
and  Hasselmann  [17]  treated  the  mode  coupling  by  the  sea-surface  roughness  and 
McDaniel  [6]  that  by  the  bottom.  Boundary  scattering  has  to  be  incorporated  in 
the  computer  code  as  a  next  step.  Also,  quantitative  measurements  on  horizontal 
sound-speed  gradients  have  to  be  made  together  with  acoustic  measurements. 

The  mode  coupling  method  yields  a  redistribution  of  intensity  in  the  vertical  plane. 
This  redistribution  is  only  dependent  on  the  sound-speed  structure,  not  on  the 
source  depth.  When  energy  is  scattered  into  the  duct  from  outside  or  vice  versa, 
the  ratio  of  the  scattered  energy  to  the  total  energy  seems  to  be  almost  frequency- 
independent.  The  total  amout  of  mode-converted  energy,  however,  is  strongly 
frequency-dependent . 

The  method  of  mode  coupling  that  has  been  used  is  a  powerful  tool  for  the  de¬ 
scription  of  random  forward  scattering.  It  is  not  dependent  on  the  mathematical 
simplifications  of  the  ray  theory,  which  prevent  direct  calculation  of  frequency- 
dependent  effects.  The  calculation  is  relatively  simple  and  fast  and  the  coupling 
calculations  are  equivalent  to  about  10%  of  th-'  mode  calculations.  In  the  same 
way  as  conventional  mode  calculations  they  provide  solutions  for  arbitrary  ranges 
without  range  dependence  penalties. 

Mode  coupling  always  affects  the  coherence  of  the  sound  field,  and  so  the  technique 
described  helps  in  defining  limits  for  source  localisation  methods  that  make  use  of 
the  modal  structure 
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RANGf  (KM) 

Fig.  1  Case  A:  comparison  of  forward  scattering  by  coupled  modes  (solid  lines)  ivit/i 
parabolic  equation  [7}  (dashed  Jines),  transmission  loss  over  range  and  depth  for  a  7  m 
deep  source  at  500  Hz 


Fig  2  Case  A  comparison  of  forward  scattering  by  coupled  modes  ( solid  lines )  with 
parabolic  equation  [7j  (dashed  lines),  transmission  loss  over  range  and  depth  for  a  7  hi 
deep  source  at  2000  Hz 
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Fig  3  Case  A:  comparison  of  /or ward  scattering  without  mode  coupling  (solid  lines)  and 
parabolic-equation  calculation  (dasbed  lines)  for  a  single  sound-speed  profile;  transmission 
Joss  over  range  and  depth  for  a  7  m  deep  source  at  oOO  Hz. 


Fig  4  Case  A  comparison  of  forward  scattering  without  mode  coupling  and  parabolic- 
equation  calculation  for  a  single  sound-speed  profile;  transmission  loss  over  range  and 
depth  for  a  7  m  deep  source  at  2000  Hz 
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5a.  Case  B  transmission  loss  at  500  He ,  calculated  (lines)  and  measured  (dots). 


Range  (km) 

5b  Case  B  transmission  loss  at  2  It  Hr.  calculated  (lines)  and  measured  (dots). 
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fig  5c  Case  B  transmission  loss  at  4  kHz,  calculated  (lines)  and  measured  (dots). 
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Fig.  6a.  Case  B.  transmission  loss  as  a  function  of  depth  with  range  with  and  without 
coupling  as  parameter;  frequency  250  Hz. 
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F’g  6b.  Case  B  transmission  Joss  as  a  function  of  depth  with  range  with  and  without 
coupling  as  parameter,  frequency  500  Hi. 
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Fif  6c.  Case  B.  transmission  loss  as  a  function  of  depth  with  range  with  and  without 
coupling  as  parameter,  frequency  J  kHw. 
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f>i  8 d  Case  B  transmission  loss  as  a  function  of  depth  with  range  with  and  without 
coupling  as  parameter,  frequency  2  kHs. 
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Fig  6e.  Case  B.  (ransmtM)on  J055 
coupling  as  parameter,  frequency  4 


as  a  function  of 
kHz 
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Fig  8a  Mode  excitation  at  range  i  rn  for  a  source  at  / '  in  depth.  frequent  y  Hz 
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F't  8b  Mode  distribution  for  15  rn  source  after  .*>(/  krn  propagation,  solid  lines  with 
mode  coupling,  dashed  lines  without  coupling 
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Fig  9e  Mode  excitation  at  1  kffi  for  deep  source  after  SO  lent  propagation,  with  coupling. 


Fig  9 f.  Mode  excitation  At  4  JtH»  for  shxUow  source  After  50  km  propAgAtion,  with  coupling 
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Fig.  10a.  Mode  spectrum  at  500  Hz  after  50  km  propagation,  with  coupling,  source 


depth  15  m. 


Fig.  10b.  Mode  spectrum  at  J  kHs  after  50  km  propagation,  with  coupling,  source  depth  50  m. 
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Fig.  10c.  Mode  spectrum  at  2  kHr  After  50  km  propagation,  with  coupling,  source  depth  15  in. 


Fig.  lOd.  Mode  spectrum  at  6.3  kHs  After  50  km  propAgAtion,  with  coupling,  source  depth  50  in. 
With  respect  to  the  computation  time  and  the  dimension  of  the  data  fields  onl.v  the  lowest  200 
modes  were  cslculsted. 
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Fig.  Jla.  Oa»*  B  transmission  toss  over 
range  a nd  depth  at  500  He.  without  mode 
coupling. 
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Fig.  lib.  Case  B:  transmission  loss  over 
range  and  depth  at  2  kHs,  without  mode 
coupling. 
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Fi*  lie.  Case  B  transmission  loss  over 
range  and  d*ptA  at  500  Hj,  with  mode 
coupling. 
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Fig  lie  C'asp  B  transmission  loss  over 
mnge  and  depth  at  2  kHz.  with  mode 
coupling,  source  beyond  the  duct  The 
steps  in  the  10 5  dB  line  are  an  artefact 
produced  by  the  limited  rpso/ution  of  levels 
of  0. 1  dB  provided  to  the  contour  program 
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Fig.  l'2a  Relative  content  of  coherent 
sound  intensity  (in  dB)  a 5  a  /unction  of 
range  and  depth;  frequency  250  Hm,  source 
depth  50  ni. 
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Fig.  12b  Relative  content  of  coherent 
sound  intensity  (in  dB)  as  a  function  of 
range  and  depth ;  frequency  250  Hs,  source 
depth  15  ni. 


5  0  10  0  15  0  20  0  25  0  JO  0  J5  0  40  0  45.0  50  0 

RANGE  (KM) 


Fig.  12c.  Relative  content  of  coherent 
sound  intensity  (in  dB)  as  a  function  of 
range  and  depth;  frequency  1  kHz,  source 
depth  50  m. 
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Fig  I2d.  Relative  content  of  coherent 
sound  intensity  (in  dB)  as  a  function  of 
range  and  depth;  frequency  1  kHz,  source 
depth  15  in. 
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Fig  J2e.  Relative  content  of  coherent 
sound  intensity  (in  dB)  as  a  function  of 
range  and  depth;  frequency  4  kHz,  source 
depth  50  m 


Fig  13.  Case  B:  sound  speed  profile 


-  <917 


I 

I 


Initial  Distribution  for  SM-199 


Ministries  of  Oefenr* 
JSPHQ  Belgium 
DND  Canada 
CHOD  Denmark 
MOD  France 
MOD  Germany 
MOD  Greece 
MOD  Italy 
MOD  Netherlands 
CHOD  Norway 
MOD  Portugal 
MOD  Spain 
MOD  Turkey 
MOD  UK 
SECDEF  US 


2 

10 

8 

8 

15 

11 

10 

12 

10 

2 

2 

5 

20 

68 


NATO  Authorities 
Defence  Planning  Committee  3 
NAMILCOM  2 

SACLANT  10 

- ' CLANTREPEUR  1 

<"  JCWE5TLANT / 

COMOCEANLANT  1 

COMSTRIKFLTANT  1 

COMIBERLANT  1 

CINCEASTLANT  1 

COMSUBACLANT  1 

COMMAIREASTLANT  1 

SACEUR  2 

CINCNORTH  ! 

CINCSOUTH  1 

COMNAVSOUTH  1 

COMSTRIKFORSOUTH  l 

COMEDCENT  ! 

COMMARAIRMED  1 

CINCHAN  3 


SCNR  Germany 
SCNR  Greece 
SCNR  Italy 
SCNR  Netherlands 
SCNR  Norway 
SCNR  Portugal 
SCNR  Turkey 
SCNR  UK 
SCNR  US 

SECGEN  Rep  SCNR 
NAMILCOM  Rep  SCNR 


National  Liaison  Officers 

NLO  Canada 

NLO  Denmark 

NLO  Germany 

NLO  Italy 

NLO  UK 

NLO  US 

NLR  to  SACLANT 
NLR  Belgium 
NLR  Canada 
NLR  Denmark 
NLR  Germany 
NLR  Greece 
NLR  Italy 
NLR  Netherlands 
NLR  Norway 
NLR  Portugal 
NLR  Turkey 
NLR  UK 


1 

1 

1 

1 

1 

1 

1 

1 

2 

1 

1 


1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

l 

1 


S C NR  for  SACLAN  TCEN 
SCNR  Belgium 
SCNR  Canada 
SCNR  Denmark 


Total  external  distribution  248 
SACIANTCEN  Library  10 

Stock 

Total  number  of  copies  280 


i 

t 

i 


40/+ 


